Bulk— Boundary Correspondence for Chiral Symmetric Quantum Walks 
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A discrete-time quantum walk (DTQW) is defined by a periodic sequence of operations on a 
quantum particle. According to the choice of the starting time of the period, different effective 
Hamiltonians can be associated to the same quantum walk. We define a DTQW to have chiral sym- 
metry (CS) if at least one of these Hamiltonians has CS. This can be ensured by using an "inversion 
symmetric" pulse sequence, which automatically gives two, inequivalent effective Hamiltonians with 
CS. We show that the sum and difference of the associated winding numbers, divided by two, gives 
the bulk topological invariants of a DTQW with CS, which control the number of and n energy 
edge states. We illustrate this bulk-boundary correspondence for the DTQW on the example of the 
"4-step quantum walk" , where tuning CS and particle-hole symmetry realizes edge states in various 
symmetry classes. 
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The realization that band insulators can have non- 
trivial topological properties which determine the low- 
energy physics at their boundary has been a rich source 
of new physics in the last decade. The general theory 
of topological insulators and superconductors [U [2] clas- 
sifies gapped Hamiltonians according to their dimension 
and their symmetries j3J. As very few real-life materials 
are topological insulators, there is a strong push to de- 
velop model systems, "artificial materials" , that simulate 
topological phases 

Discrete-time quantum walks (DTQW) [SHE] offer a 
versatile way of simulating topological insulators [SHU]- 
Beyond that, they also have unique topological fea- 
tures, with no counterpart in standard solid-state se- 
tups. In DTQWs with particle-hole symmetry (PHS), 
edge states, "Majorana modes" can have two protected 
quasienergies: e = or ir (time is measured in units of 
the timestep and h = 1). Generalizing the results for 
periodically driven systems [12] , a pair of topological in- 
variants (Z 2 x Z 2 ) has been defined, which control the 
appearance of these topologically protected bound states 
[13] , Both and n energy Majorana edge states have 
been experimentally observed in an optical realization of 
a quantum walk |14] . 

As for chiral symmetry (CS), the situtation of 
DTQW's is much less clear. Even for the simplest 1- 
dimensional DTQW, it is disputed whether it even has 
CS [5] or not [12] • Although it is expected that CS should 
imply aZxZ bulk topological invariant, this has not 
yet been found for DTQW's. As opposed to the case of 
PHS, there is also not much to draw on from periodically 
driven systems. What DTQW's have CS? How can the 
bulk "winding number" be expressed for DTQW's with 
CS? These are the problems we tackle in this paper. 

A DTQW concerns the dynamics of a particle, 
"walker", whose wavefunction is given by a vector, \ty) = 



Ex^Es^i.i*^*)! 1 . 5 )- Here ; x = f,-.., A is the 
discrete position, and s = ±1 indexes the two orthogo- 
nal internal states of the walker, the "coin eigenstates" , 
which we also refer to as "spin". The dynamics, instead 
of given by a time-independent Hamiltonian, is realized 
using a periodic sequence of alternating "step" and "coin 
rotation" operations. 

The step operations are translations of the particle by 
one lattice site depending on the value of the "coin" , the 
z-component of its spin. These are described by unitary 
operators S s , where s is either + or — , and 

N 

S±=X;(|s±1><s|®I±1><±1| + \x){x\®\Tl)(Tl\). (1) 

x=l 

For simplicity, we take periodic boundary conditions. 

Between each two steps, a site-dependent local "coin 
rotation" Rj on the walker's internal state is performed. 
We consider 

N 

R j = J2\z)(x\®R(x j (x),o j (x)); (2) 

x=l 

, n\ _ f cos8 — ismOs'mx — sin cos x \ / q \ 
^ ' \ sin 6 cos x cos + i sin 8 sin \) 

= exp[-i0(cos(x)ci, + sm(x)<J z )} (4) 

This allows breaking PHS via the angle x 0- Details 
of how the local operations Rj are performed do not in- 
fluence the DTQW, all the information about them is 
summarized in the corresponding unitaries Rj. 

One period of the DTQW is defined by |W(t + 1)> = 
U \^>(t)), for t € Z. Here the unitary timestep (Floquet) 
operator is composed of 2M successive pulses, 

Uo = SmRmSm-iRm-i ■ ■ ■ S\R\. (5) 

A period has to include an equal number of 5+ and 
5_ pulses, otherwise timestep operator has quasienergy 
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winding [T5], and cannot have gaps. Thus, M is even. 
We take each pulse to have a duration 1/2M, with- 
out losing generality. We can also shift the start- 
ing time of the period by T, giving the dynamics as 
\^(t + 1+T)) = U T \^(t + T)), for any t G Z. The start- 
ing time T has to be during a rotation, restricting it to 
T = (I - 1)/M + y/(2M), for 1 < I < M, and < y < 1. 
The shifted Floquet operator reads 

Ut = R\i + iSm+i-\Rm+i-i ■ ■ ■ Si+iRi+iSiR] v , (6) 

where R 1 - = Y^ x \ x )( x \ ® R(Xj{ x )> V 9j{ x ))- This is a uni- 
tary transform of Uq, we refer to the transformation con- 
necting the Ut to Uq as a "gauge" tranformation. Wc 
don't consider splitting translation operators S. In many 
realizations of DTQWs, e.g. the ones on atoms trapped 
in optical lattices or on photons on an optical table, per- 
forming only part of a shift operation would leave the 
walker between sites, and its description would necessi- 
tate an increased Hilbert space. 

A DTQW can be seen as a stroboscopic simulator of 
an effective Hamiltonian _ff c ff.T, which is associated to 
the Floquet operator by 



(7) 



The effective Hamiltonian is uniquely defined by this 
equation if we restrict its eigenvalues, the quasienergies, 
to an "energy Brillouin zone", — tt < e < tt. This is 
completely analogous to the restriction of the quasimo- 
mentum to the first Brillouin zone. 

Previous work on CS in DTQWs has focussed on 
whether -ff e ff,o has CS, and identifying the associated 
topological invariant. We realize that it is important to 
widen the scope: a DTQW has CS, if there is a gauge 
where its effective Hamiltonian has CS: If a time T and 
a unitary operator T acting on the coin space can be 
found, with T 2 = 1, that TU T T = U^ 1 . This unitary 
time-reversal symmetry is equivalent to CS of the effec- 
tive Hamiltonian i? ff,T- 

We now proceed to derive the bulk-boundary corre- 
spondence for DTQWs with CS, i.e., a formula for the 
Z x Z bulk topological invariant controlling the appear- 
ance of and tt energy edge states. To do this, we first 
show that DTQWs with CS have not only one, but two 
inequivalent chiral symmetric gauges. 

Based on the definition of CS above, a DTQW has 
CS represented by T — a x , if the sequence of operations 
defining the walk has an "inversion point" . By this, we 
mean that there is an I, with which for every j: 



Ri-j — Ri+j\ 

Si-j = S+ O Sl+j+l = S-. 



(8) 
(9) 



As we show below, this not only ensures that there is a 
gauge where the effective Hamiltonian has CS, but two 
inequivalent such gauges. 



To prove that "inversion symmetry" of the pulse se- 
quence gives two CS gauges, we begin by noticing that 



whereby also o~ x S+<j x 



(10) 

Furthermore, since the 
local unitaries are rotations Rj about axes that have no 
x-component, 

°~ x R(0,x)&i 



Rie^r 1 = R(-e,x)- 



(ii) 



Take the sequences of M operations just after and just 
before the middle of the "inversion point" , 



F = R 1 1 +~ m / 2 Si-i+m/2Ri-i+m/2 ■ ■ ■ Ri+iSiRi'^; (12) 

G = R.y 2 Si-iRi-i . . . Ri +1 -m /2S1-M i%R\-mi%' (13) 
These give us two Floquet operators for the walk: 

U' = FG; U" = GF. (14) 



1/2. 



Using relations ( 10 ) and ( 11 ), we have that time reversal 



can be done during a period, TFTG = 1, whereby G — 
rF _1 r. From this it is straightforward to show that both 
U' and U" are chiral symmetric. 

CS allows a definition of sublattices, via the projection 
operators U A = (1 + T)/2, U B = (1 - T)/2. Eigenstates 
of H' ee with quasienergy e 7^ 0, tt can be chosen to have 
equal support on both sublattices. Stationary states with 
quasienergies or tt, however, can be chosen to be on a 
single sublattice in the chiral gauge U' , i.e., their wave- 
functions in this gauge are eigenstates of T. 

For the bulk-boundary correspondence, we consider an 
inhomogeneous DTQW with CS, consisting of a transla- 
tionally invariant "L" bulk at 1 < a; < d and an "i?" 
bulk at d <C x <C N. There are (smooth or sharp) bound- 
aries between the two bulks around x ~ d and x w 1. In 
the chiral symmetric gauge U', the two bulks have ef- 
fective Hamiltonians H' cS L and H' eSR . We assume both 
bulk Hamiltonians have gaps around e = and e = tt. 
Therefore, stationary states with quasienergies e = or 
tt must have wavefunctions confined to the edges, expo- 
nentially decaying towards the bulks. The edge around 
x s» d hosts m' A (m' B ) edge states on sublattice A (B), 
given by 



*A = m A,a + m 'A,n\ m 'B = 



l B,0 



L B,TV1 



(15) 



where the second index stands for the energy. We are 
looking for the topological invariants of the bulk parts 
of the walk, vx,ei with X — L, R and e = 0, tt, whose 
differences give us the number of topologically protected 
edge states separately for each energy, 



l B,e- 



(16) 



The first step towards the topological invariants is the 
standard winding number v' x [5] associated to the bulk 
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effective Hamiltonian H' cS x , where X = L or X = R 
in gauge U' . This is obtained from the bulk, transla- 
tional invariant part of the Floquet operator, diagonal 
in momentum space: U' x = ^2 k \k)(k\ <S> U' x (k), and 
U' x (k) = e - iH ° a ' x< - k \ Instead of the effective Hamil- 
tonian, it is convenient to calculate with H' x (k) — 
sva{H C K : x{k)}- This has the same CS and the same wind- 
ing number as H e ^x{k) (and as its flattened version 
Q = sgn{H c ^ x (k)] 0), but can be obtained much more 
efficiently via H' x (k) = [U' x (k)^ -U' x {k)]/(2i). In a basis 
where V = diag(l, . . . , 1, —1, . . . , —1) is a diagonal matrix 
with an equal number of +1 and —1 elements, the matrix 
of H' x (k) is block off diagonal because of CS. We name 
its upper right block h' x [k). The winding number v' x 
reads 



We compare this with Eq. (16), and simply read off the 



'x 



1 



dk — log det h'x(k). 
ak 



(17) 



The winding number v' x is related to the difference 
of the bulk polarizations on the two sublattices in bulk 
X [IB]. Therefore, it cannot differentiate between and 
7r energy edge states, and can only be used to obtain 
the sum of all topologically protected edge states around 
x « d: 



'L - - "M.,0 "I" m A,7T - m B,0 - m B,n- ( 18 ) 

However, there is the other CS gauge, U", where we have 



'R 



m A,0 + m A, 



l B,0 



' B. 



(19) 



We need to combine the information from the two CS 
gauges to obtain the topological invariants. 

We can obtain a simple connection between the two CS 
gauges by considering an edge state. In the gauge of U' = 
FG, the edge state has a wavefunction 'J, entirely on 
sublattice A (or B), i.e., = (-1) 9 *, with g = (or 1). 
In other gauges, where Ut has no CS, the energy of the 
edge state has to remain the same, but its wavefunction 
can extend over both sublattices. In the other CS gauge 
U" = GF, however, its wavefunction, <f> = G*ff, again has 
to to be entirely on a single sublattice. This can be A (or 
B), whereby T$ = (— 1)^ with / = (or 1). Consider 

GF& = GrG _1 r$ = gtg- 1 (-i) / * = Gr(-l)^ = 
(-l)B+fG$ = {-l) 9+f $. This shows that ( 7r) energy 
edge states are on the same (opposite) sublattice in the 
two CS gauges. This can be summarized as 



l A - '"71,0 1" '">B,m m B ~ m B,0 + m A,7T- ( 20 ) 

To obtain the number of protected edge states at zero 



and 7r energies separately, we simply substitute Eqs. ( 20 ) 



into Eqs. (|18|), (j 1 5[) and (|19), and rearrange to obtain 

■ - v" v' -v" 

_ (22 ) 



m' A - m! 



bulk topological invariants (vq, v w ) as 



(^o,^tt) = 



V + V v — V 



(23) 



This, the bulk-edge correspondence for DTQWs with 
CS, is the main result of the paper. 

Having derived a general formula for the topologi- 
cal invariant of ID DTQWs with CS, we now discuss 
an example where the differences between CS and PHS 
come into play. To arrive to the example, first con- 
sider the "split-step walk" of Kitagawa et al. [9] , given by 
G =5+^(0,^)5-^(0,01). There, both PHS and CS are 
present, and we find that the v e are one-to-one functions 
of the invariants Q £ induced by PHS [13 : v E = 1/2 — Q £ , 
for both e = 0, ir. (An interesting special case is the 
simple quantum walk, obtained by setting = 0.) We 
can break PHS by using nonzero angles \- To be able 
to break CS, we can consider a longer period of pulses, a 
"4-step DTQW" , given by 



Uq — S-^R^S^-R^S— R2S— R±. 



(24) 



This walk has no CS if R 2 # R4, but has CS if R 2 = R 4 , 
with F = rI ,2 S-R 2 S-R\' 2 and G = R\ /2 S+R 4 S+rI /2 . 
The 4-step walk also has the advantage that the effective 
Hamiltonian will have longer range hoppings, and thus we 
can expect higher values of the winding numbers. This 
is entirely analogous to adding a 3rd nearest neighbor 
hopping term to the SSH model. 

The topological invariants in a section of the phase 
space of the 4-step DTQW with both CS and PHS (Car- 
tan class BDI [3]) are shown in Fig. [I] Here we set all 
Xj = to ensure PHS, 9 2 = #4 to ensure CS, and set 
#i=0 for simplicity. We restrict 9 2 to — 7r/2 < 9 2 < n/2 
since adding 77 to both 9 2 and #4 just brings two factors of 
— 1 that cancel out in both chiral gauges. Generic values 
of 9 2 = #4 and $3 give effective Hamiltonians with gaps 
around both e = and e — n. Examples for these are the 
points C(9 2 = tt/20,03 = tt/4), D [9 2 = 0, 9 3 = -tt/4), 
and E(9 2 = n/4,9 3 = tt/4). 

To see the effects of breaking the symmetries on edge 
states, we consider two inhomogeneous systems, consist- 
ing of two domains of 40 sites each, with sharp boundaries 
in between. The inhomogeneous rotations read 



40 



SO 



Rj = \ x )( x \ ® Rj,x + I s ) ( x \ ® ( 25 ) 

x=l x=41 

where X = D or X = E, and G refer to the parameter 
sets of defined in the previous paragraph and indicated 
in Fig. [IJ 

We break PHS (realizing Cartan class AM) in a con- 
trolled way by introducing a nonzero X3 i n the bulk 
< x < 41. As long as the bulk gaps are still open, 
breaking PHS does not change the edge state energies, 




FIG. 1: (Color online) Parameter space of the 4-step DTQW 
with PHS ensured by \j = 0, CS ensured by 82 = 84, and 81 
set to 0. The DTQW has effective Hamiltonians with gaps 
around both e = and e — n, except at the gapless points 
where gaps close at e — (solid lines) or e = tt (dashed 
lines). Single lines indicate that the gap closes at a single 
k, at either k — or k — ir. Double lines indicate double 
gap closings, at k — ±7r/2. For each gapped domain, the 
corresponding pair of winding numbers {u',i>"} as well as 
the pair of topological invariants (yo, z-v), cf. Eq. ( |23[ ), are 
shown. Letters "C" , "D" and "E" indicate sets of parameters 
used for the inhomogeneous quantum walk, with rotation as 
in Eq. p5l. 



as shown in Fig. [2]A1), Bl). The edge state energies are 
still protected by CS, and can only move from their origi- 
nal values if the bulk gap closes (at X3 = 7r /2 for the D-C 
boundary). We break CS (realizing Cartan class D), by 
changing 62 — O4 in the "L" bulk. A pair of edge states on 
the same edge at the same energy can now break apart, 
becoming PHS partners of each other. This can be seen 
in Fig. [2] A2) at both and tt energy. However, a single 
edge state, as the one between bulks B and C, is still 
protected by PHS when CS is broken, as seen in Fig. [2] 
B2). Finally, to check that no extra hidden symmetries 
remain, we break both CS and PHS (realizing Cartan 
class A). In that case the edge state energies are not pro- 
tected anymore, cf Fig. [2] A3), B3). This shows that our 
description of the relevant symmetries of the DTQW was 
indeed exhaustive. 

To summarize, we gave a definition of CS for DTQWs, 
and derived the corresponding bulk topological invari- 
ants, using the fact that the walk is defined by a sequence 
of operations, rather than just by its unitary timestep 
operator. The approach presented here based on finding 
the "symmetric gauges" should generalize to periodically 
driven quantum systems [T51 [PTH2TJ] . In such setups, PHS 
has been shown to lead to 2 types of "Floquet Majorana 
fermions" [12], which should have clear signals in trans- 
port |21j and can also be useful for quantum information 
processing [55]. Theoretical proposals have already seen 
several such states at a single edge if the driving also en- 
sures CS [33]. The bulk topological invariant controlling 




3it/2 2jt kIS 71 371/2 2% 

x 3 =(e 2 -e 4 )/4 



FIG. 2: (Color online) Spectra of an inhomogeneous "4-step 
walk" on N = 60 sites as defined by Eqs. pi|,p5]|, with 
two domains: D and C (top row), and E and C (bottom 
row). Left panels: We break PHS via X3 for x < 41 on Al) 
and Bl). As long as the bulk gaps are still open, the edge 
state energies do not change as they are still protected by 
CS. Middle panels, A2) and B2): we break CS by setting 
8 2 -> 8 2 + A8, 84 -> 84 - A8 for n < 31. This lifts the 
degeneracy of edge states on the same edge pairwise. At the 
interface between E and C, the unpaired edge state remains 
(B2). Right panels, A3) and B3): as both PHS and CS are 
broken, no topologically protected edge states remain. In B3) , 
the unpaired edge states at both edges are displaced in energy. 



the number of these edge states is as yet unknown, but 
it could be derived using the approach of this paper. 

Although 2-dimensional DTQWs have already been re- 
alized in experiments [24) , their topological invariants are 
largely unexplored. In 2 dimensions, edge states can exist 
in the absence of symmetries; the related bulk-boundary 
correspondence for periodically drive systems has only 
recently been found |25j . The approach of identifying 
the "symmetric gauges" should be a key idea for the de- 
scription of other symmetry classes for both periodically 
driven quantum systems and DTQWs. 
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